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Biochemistry and mechanics are closely coupled in cell adhesion. At sites of cell- 
matrix adhesion, mechanical force triggers signaling through the Rho-pathway, which 
leads to structural reinforcement and increased contractility in the actin cytoskeleton. 
The resulting force acts back to the sites of adhesion, resulting in a positive feedback 
loop for mature adhesion. Here we model this biochemical-mechanical feedback loop 
for the special case when the actin cytoskeleton is organized in stress fibers, which 
are contractile bundles of actin filaments. Activation of myosin II molecular motors 
through the Rho-pathway is described by a system of reaction-diffusion equations, 
which are coupled into a viscoelastic model for a contractile actin bundle. We find 
strong spatial gradients in the activation of contractility and in the corresponding 
deformation pattern of the stress fiber, in good agreement with experimental findings. 



I. INTRODUCTION 

Adhesion of biological cells to each other and the extracellular matrix is one of the hall- 
marks of multicellular organisms and a very active area of research in cell biology. Investi- 
gating cell adhesion is crucial to understand physiological processes like tissue development 
and maintenance, but also disease-related processes like growth and migration of cancer cells 
. In general, biological cells have a limited repertoire of possible behaviours, including 
spreading, adhesion, migration, division or death, but very sophisticated ways of controling 
the switching between these different processes. Traditionally, investigation of this cellular 
decision making has focused strongly on the biochemical aspects, including detailed models 
for signal transduction {^,4]. During recent years, it has become increasingly clear that in 
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adhesion-related processes, cellular behaviour is not only controlled by biochemical cues, but 
also involves many physical determinants like the structural organization of the extracellu- 
lar matrix and the cytoskeleton or force generation through molecular motors aaaa. 
For example, it has been shown that the stiffness of the extracellular environment deter- 
mines migration of tissue cells and different at ion of stem cells 10 1. In particular, these 
celluar responses have been found to depend on the ability of the cells to contract their 
environment with actomyosin contractility and to convert this mechanical process into a 
biochemical signal. Although these processes are essential for such important situations like 
tissue functioning or cancer cell migration, theoretical models describing the coupling be- 
tween biochemistry and mechanics in cell adhesion are still rare, albeit essential for a future 
systematic understanding of how multicellular organisms organize themselves. 

Cell adhesion is closely related to the actin cytoskeleton, whose organization is central 
in determining the structural properties of cells. In cell culture with stiff substrates, the 
actin cytoskeleton tends to organize in stress fibers, which are bundles of actin filaments 
tensed by myosin II molecular motors Uj. Stress fibers usually end in focal adhesions, 
which are integrin-based adhesion contacts which can grow to a lateral size of several microns 
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13| . On their cytoplasmic side, focal adhesions recruit more than 90 components (mostly 
proteins) which physically reside in the adhesion structure [l^. In 1992, Ridley and Hall 
published a landmark paper demonstrating that the assembly of stress fibers and focal 
adhesions is regulated by a small GTPase called Rho [isl. Rho has many isoforms, but the 
one mainly associated with focal adhesions is RhoA, which for simplicity in the following 
we refered to as Rho. In a companion paper of the same year, Ridley and Hall showed 
together with coworkers that another small GTPase called Rac stimulates the formation 



of lamellipodia as they appear in cell migration pj|. The main isoform associated with 
focal adhesions is Racl which for simplicity in the following we will refer to as Rac. While 
Rho mainly acts through activation of actomyosin contractility, the main effect of Rac is 
activation of actin polymerization, in particular activation of the actin nucleation factor 
Arp2/3. It has been reported later that activation of Rac downregulates Rho, leading to 
disassembly of stress fibers and focal adhesions [l^]. In many situations, Rho and Rac can 
)e regarded as antagonists, switching the cytoskeleton between different structural states 
isl . They are part of a larger family of small GTPases, called the Rho-family, which for 
example also includes Cdc42, which stimulates the formation of filopodia and maintains 



cell polarity [19|]. Apart from regulation of the actin cytoskeleton, the small GTPases from 
the Rho-family have many other functions in the cell, for example in cell cycle control and 
differentiation. 

Although the small GTPases from the Rho-family are simple molecular switches, they 
are regulated by many different factors. In general, GTPases are upregulated by guanine 
nucleotide exchange factors (GEFs), which convert the inactive Rho-GDP form to the active 
Rho-GTP form by exchanging GDP for GTP. They are downregulated by GTPase-activating 
proteins (GAPs), which stimulate Rho-GTPase activity, thus leading to GTP-hydrolysis 
and transforming the active Rho-GTP to the inactive Rho-GDP. For the 20 members of 
the Rho-family, 60 different GEFs and 70 different GAPs as well as more than 60 different 



downstream targets have been identified 20(]. At the current stage of affairs, there is no 



way how this complex network can be modeled in full detail. However, certain parts of this 
network have been well characterized by biochemical assays, in particular different parts of 
the Rho-mediated signal transduction pathway leading from focal adhesions to actomyosin 
contractility. 

Here we focus on the role of Rho as stabilizing factor for mature adhesion. During recent 
years, it has been shown that Rho is the central component of a biochemical-mechanical 
feedback loop which regulates mature adhesion. In detail, it has been shown that application 



of force on focal adhesions triggers their growth in a Rho-dependent manner 
in 
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21| (reviewed 



181]). Two main downstream targets of Rho leading to stress fiber formation have 
been identified. The formin mDia leads to actin polymerization, while the Rho-associated 
kinase ROCK leads to phosphorylation of myosin light chain and thus to increased motor 
activity. Together these effects lead to formation of and contractility in stress fibers and 
therefore to increased force levels at focal adhesions. In this way, a positive feedback loop is 
closed for upregulation of mature adhesion characterized by focal adhesions and stress fibers. 
This biochemical-mechanical feedback loop is schematically depicted in Fig. [TJ An essential 
part of this feedback loop is the growth of focal adhesion under force, which recently has 



been the subject of different modelling approaches [22, l23|, l2J, |25|, l26|, |27| (reviewed in (28| ) . 
However, these models have focused mainly on the mechanical and thermodynamic aspects of 
the growth process, neglecting the interaction of mechanics and biochemical signaling. The 
positive feedback loop between contractility and growth of adhesions has been modelled 
before in the framework of kinetic equations, but without addressing the details of force 
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generation and its regulation by signaling pathways 29| . Similar kinetic equations have been 
used to model the antagonistic roles of Rho and Rac in cell adhesion, but again without 
addressing the details of force generation and regulation 30[. Recently force generation 
has been addressed in more detail in a model for whole cell contractility and stress fiber 



formation 
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32l | . However, no details of the signalling pathway have been modelled except 



for an unspecified activation signal. 

Because the actin cytoskeleton is very dynamic and interacts with many different molec- 
ular factors, including actin-associated proteins and molecular motors, it is very difficult to 
model its mechanical properties in a general way. However, modelling becomes feasible if 
one focuses on one of the well-characterized states of the actin cytoskeleton, for example the 
lamellipodium or stress fibers. Because here we are mostly interested in mature cell adhe- 
sion in culture, we will focus on the latter case. Modelling stress fibers can be approached 
from different perspectives. An obvious starting point are their common characteristics with 
muscle fibers, which is a linear sequence of sarcomeres, each containing around 300 myosin 
II molecular motors working collectively together as they slide the actin filaments relatively 
to each other. This field has been pioneered by the Huxley-model which later has 



been modified in many regards, e.g. in regard to filament extensibility [3J] or by a detailed 
modelling of the myosin II hydrolysis cycle 35|]. In contrast to muscle fibers, stress fibers 
are more disordered and a complete description therefore requires a model for their assem- 
bly process from polar filaments interacting through molecular motors. Such a description 
has been achieved in the framework of a phenomenological theory which however does not 
model the details of the underlying motor activity 36|, l37|]. This theory does predict differ- 
ent dynamical states of the system, including a stationary state of isometric contraction as 
observed in stress fibers. 

Although being less ordered than muscle on the level of electron microscopy, stress fibers 
do exhibit a periodic organization. Fig. [2] shows experimental data for fibroblast adhesion to 



a stiff substrate 38|]. The image of a whole cell shown in Fig. [2JA. reveals a banding pattern 
in the stress fibers. The green regions correspond to the actin crosslinker a-actinin while the 
red regions correspond to myosin II molecular motors. Non-muscle myosin II is known to 
assemble into bipolar minifilaments consisting of 10-30 myosins, as depicted in the cartoon 
of a stress fiber in Fig. [TJ In contrast to muscle, the striation pattern of stress fibers shows 
considerable variability along the length of a stress fiber. In particular, it has been observed 
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that upon stimulation of contraction with the drug calycuhn A, only the sarcomeres in the 
periphery (close to the focal adhesions, Fig. [2l3) shorten, while those in the center (close to 
the cell body. Fig. [2lC) elongate. The exact time course of the striation width is shown in 
Fig. [2p and E for control and stimulation experiments, respectively. This demonstrates that 



activation of contractility leads to strong spatial gradients in the striation patterm 
not sufficient to model only one sarcomeric unit as it is typically done for muscle 



33 



lus it is 
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35|. 



Rather at least a one-dimensional chain of such sarcomeres has to be considered. 

Here we introduce a new one-dimensional model for stress fibers which essentially models 
them as linear sequence of viscoelastic Kelvin- Voigt bodies, whose stationary state is deter- 
mined by the elastic part. The action of the molecular motors inside each sarcomeric unit 
are included on the level of a linearized force-velocity relationship. The signaling pathway 
is modeled as a system of react ion- diffusion equations. For this study, we have conducted 
an extensive survey of the relevant literature and have collected the measured rate and 
diffusion constants in such a way that they now can be used for mathematical modelling. 
The coupling between the biochemical signaling pathway and the mechanical stress fiber 
model proceeds by introducing a spatially varying fraction of active molecular motors. The 
local activation level is thereby determined by the outcome of the signaling pathway. A 
continuum limit of the mechanical model for many sarcomeric units in series results in a 
partial differential equation with mixed derivatives. The whole system of reaction-diffusion 
equations for the signal transduction and the partial differential equation for the mechanical 
part can be solved simultaneously. By feeding the force resulting from the mechanical model 
back into the activation of the signaling pathway, we obtain for the first time a model for 
the closed biochemical-mechanical feedback cycle described above. As a first application of 
our model, here we show that it predicts heterogenous contraction of stress fibers, in good 
agreement with experiments. In general, our work shows how models for the coupling of 
biochemistry and mechanics can be deviced in a meaningful way. 

The paper is organized as follows. In Sec. [Ill we start with our model for the part of the 
Rho-pathway which is relevant for our purposes. We then introduce our mechanical model 
for one sarcomeric unit of the stress fibers in Sec. IIIII and its continuum limit for a whole 
stress fiber in Sec. IIVI In Sec. IVlour model predictions for inhomogeneous contraction upon 
activation of contractility by calyculin are compared to the experimental results. Finally we 
conclude in Sec. IVD 
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II. BIOCHEMICAL MODEL FOR THE RHO-PATHWAY 

In this study we concentrate on stress fibers and their regulation by the Rho-pathway. 
In Fig. [U we introduce a coordinate system for our one-dimensional model: the stress fiber 
extends along the positive x-direction and the endpoints at a; = and x = L correspond 
to two focal adhesions. Because the two focal adhesions are treated as equivalent, our 
model has inflection symmetry around x = L/2. In Fig. [3l we schematically depict the 
biochemical part of our model in the spatial context of the focal adhesion at x = (by 
symmetry, the same description applies to the one at x = L). Three compartments have to 
be considered: the focal adhesion, the cytoplasm and the stress fiber. In our model, each 
of these compartments corresponds to one or two important biochemical components. The 
reaction pathway is a linear sequence of activating or inihibitory enzyme reactions initiated 
at focal adhesion, transmitted through the cytoplasm by diffusion and resulting in spatially 
dependent myosin activation in the stress fiber. In the following we discuss each reaction 
step in detail and show how these processes are translated into reaction-diffusion equations. 
The abbrevations used for the biochemical components are compiled in Tab. [H together with 
a short description of their functions. The model equations are summarized in Tab. [Tll 

Our modelling of the biochemical signaling pathway starts with the activation of Rho at 
focal adhesions. The Rho-protein has a lipophilic end that serves to anchor it to lipid mem- 
branes [3^. Complexation with guanine nucleotide dissociation inhibitors (GDIs) shields 
the hydrophobic parts of the Rho-protein and make it inactive as well as soluble in the 
cytoplasm {40]. It is expected that Rho is released from these complexes at focal adhesions. 
More importantly, focal adhesions are known to recruit different Rho-GEFs, thus activating 
Rho at the focal adhesions. The active Rho is then able to bind Rho-associated kinase 
(ROCK) and thereby activate its kinase activity 4l|. Since the active ROCK is bound to 
Rho-GTP, we assume in our model that these components are not diffusible but are local- 
ized to the focal adhesions. For the same reason, we neglect the direct interaction of ROCK 
and myosin which has been reported to occur in vitro. Active ROCK phosphorylates the 
diffusible myosin light chain phosphatase (MLCP) at its myosin-binding subunit (MBS). 
MLCP and its antagonistic partner myosin light chain kinase (MLCK) are the main reg- 
ulators of myosin contractility in the context we are interested in. Both enzymes interact 
with the regulatory myosin light chain subunit (MLC) of the myosins. Depending on the 
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respective activities, MLC gets either phosphorylated or dephosphorylated. MLC, in turn, 
controls the myosin binding to actin filaments. Only if MLC is phosphorylated myosin is 
able to bind actin filaments and perform its ATPase cycle that converts chemical energy into 
mechanical work, e.g. contraction of stress fibers 421]. By phosphorylating MLCP, ROCK 
effectively enhances the phosphorylation level of MLC. In this way, Rho-activation can lead 
to an increase in myosin contractility in the stress fiber. 

Our model for the biochemical react ion- diffusion system assumes that each enzyme stim- 



ulation follows Michaelis-Menten kinetics 



43|. In Michaelis-Menten kinetics, production 



first increases linearly with educt concentration and then saturates at a maximal production 
velocity if educt concentration exceeds the value set by the Michaelis-Menten constant. The 
precise molecular process corresponding to the conversion of force into a biochemical signal 
at focal adhesions has not been identified yet. However, it is expected that mechanical forces 
exerted onto focal adhesions eventually initiate the loading of Rho-GTP leading to ROCK 
activation. For lack of information we therefore lump the focal adhesion associated pro- 
cesses into one equation that effectively describes the conversion of ROCK into its activated 
form (presumably complexed with Rho-GTP). The mechanical force Fb that stimulates the 
activation is treated as enzyme in the framework of Michaelis-Menten kinetics: 

dROCKjt) _ riFb{t){ROCKtot - ROCKjt)) V_iROCK{t) 

dt ~ Ki + {ROCKtot - ROC Kit)) " K_i + ROCK{t) ^' 
The variable ROCK denotes the activated form of ROCK and we assume that the overall 
concentration of ROCK is constant at ROCKtot- The force exerted by the stress fiber 
onto the focal adhesion, Ff,{t), stimulates the conversion of ROCK into its activated form 
with maximum velocity riFh{t) and Michaelis-Menten constant Ki. The parameter ri is 
equivalent to a rate constant but relates mechanical force to a chemical reaction. For this 
reason the units of ri are given as [nM/snN]. The force Fh{t) will depend on the stress 
fiber deformation. The second term accounts for the degradation of activated ROCK to its 
inactive form, with maximum velocity V-i and Michaelis-Menten constant K^i. Since we 
expect ROCK in its active form to be associated with focal adhesions, we omit diffusive 
contributions to this equation. 

One main effector of ROCK is MLCP which we regard as a diffusible compound leading 
to a reaction-diffusion equation: 

dMLCP{x,t) _ V_2MLCF-P{x,t) d'^MLCP{x,t) 
di ~ K_2 + MLCP-P{x,t) ^ dx^ 
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Here, the variables MLCP-P and MLCP denote the phosphorylated and unphosphorylated 
form of myosin hght chain phosphatase, respectively. The first term accounts for the de- 
phosphorylation of MLCP-P with maximum velocity V-2 and Michaelis-Menten constant 
K^2- The second term allows for the diffusion of the phosphatase with diffusion constant 
D. The phosphorylation level of MLCP is also regulated by the active form of ROCK which 
catalyzes the reverse reaction, that is the conversion of the phosphatase into its phosphory- 
lated form. However, ROCK is only active in the vicinity of focal adhesions located at each 
end of the stress fiber. Therefore this source term can be incorporated into the boundary 
conditions for Eq. ([2]), in the sense that the diffusive flux into the boundary has to balance 
the conversion into its inactive form: 

dMLCPjx = 0,t) _ R2R0CK{t)MLCP{x = 0,t) 

dx ~ K2 + MLCP{x = 0,t) ^ ^ 

The same relation, but with inverted sign is valid at the other end at x = L, compare the 
overview in Tab. [Ill The reaction is again modelled with Michaelis-Menten kinetics, where 
i?2 = '^2Vb is the product of a rate constant r2 with an effective volume vt, of the focal 
adhesion in which the reaction takes place. K2 is the usual Michaelis-Menten constant. For 
the phosphorylated form of the phosphatase (MLCP-P) equations similar to Eq. ([2]) and 
Eq. (IH]) are valid, but with inverted signs of the source terms and a diffusion constant which 
in principle can be different, see Tab. [TTl 

MLCP together with MLCK regulate the phosphorylation level of MLC. Since myosin in 
stress fibers form mini- filaments which are bound to actin filaments, we neglect diffusion of 
this compound, leading to the rate equation for the phosphorylated fraction n of MLC: 

dn{x,t) _ Vsil -n{x,t)) r.^MLCP{x,i)n{x,t) 
di ~ K^ + {1 -n{x,t)) K_s I + n{x,t) 

By allowing only the ratio of the phosphorylated fraction to vary, we assume that the overall 
amount of myosin in the stress fibers is fixed. MLC is phosphorylated by MLCK with a 
maximum velocity V3 = r^MLCK and respective Michaelis-Menten constant K3. Here we 
assume that the concentration of MLCK is constant within the cell. The kinase is antago- 
nized by MLCP that dephosphorylates MLC with a rate constant r_3 and Michaelis-Menten 
constant A'_3. The factor I is an inhibition parameter defined below. Since MLCP has 
spatial dependent source terms and is diffusible, the inhibition of MLC by the phosphatase 
will vary in space. 
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To complete the biochemical modelling we have to specify how the induction of calyculin 
is treated in our model. Calyculin is an inhibitor of MLCP and thereby enhances the phos- 
phorylation level of MLC. We model the interaction of calyculin with its target MLCP as a 



competitive inhibition leading to the additional factor / in the last term of Eq. (j4]) [4J] . In 
the presence of calyculin / > 1 (in absence of the drug: / = 1) which effectively increases 
the Michaelis-Menten constant and thus decreases the rate of MLC dephosphorylation. 
Hence more myosin motors will be activated and cell contractility is stimulated. The in- 
duction of the drug is then modelled by switching instantaneously the inhibition parameter 
from J = 1 to J = 3. Thereby we omit the time delay caused by the internalization of the 
drug. 

The used parameter values for the reaction-diffusion system are based on an extensive 
survey of the literature and are summarized in Tab. Illli If a range of values is reported 
in the literature, we chose an intermediate value for this parameter. If no value could be 
found in the literature, we made reasonable assumptions based on similar parameters in 
other systems. No attempt was made to fit the parameters to some target function. We 
first analyze the properties of this reaction diffusion system assuming that the boundary 
force exerted onto the focal adhesions is held at a constant level. This would be the case 
if the myosin forces were in a stationary state and not regulated by the biochemical signals 
emerging from focal adhesions. We impose artificial initial conditions that all components 
{ROCK, MLCP-P and n) are at zero activation level but set the boundary force to 5 nN, a 
typical force observed for fibroblast 45|]. This mechanical stimulation triggers the accumu- 
lation of active ROCK at focal adhesions, see Eq. ([T]), creating a sink for the active form of 
MLCP. Thus in those boundary regions MLCK dominates and increases MLC phosphory- 
lation. Closer to the center of the cell MLC rather remains in its unphosphorylated form. 
The width of the interfacial region of intermediate MLC-phosphorylation level is mainly 
determined by the diffusiveness of the phosphatase. The faster the diffusion, the wider the 
intermediate region. On a typical timescale of a few minutes all components equilibrate to 
their steady state concentration profile, where MLC is highly phosphorylated at the bound- 
aries, but is poorly activated at the center of the stress fiber. In Fig. H] we show the typical 
equilibration of the phosphorylated fraction of MLC, t), as obtained from a solution of 
the full system of biochemical react ion- diffusion equations. Below we will argue that this 
phosphorylation profile of MLC implies a spatially varying myosin motor activation leading 
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to an inhomogeneous stress fiber contraction. 

III. SARCOMERIC UNIT OF STRESS FIBER MODEL 

A minimal model for stress fibers has to take into account not only the viscoelastic but 
also the contractile properties of the fiber due to myosin motor activity. For the mechanical 
resp onse of a sarcomeric unit we take the usual Kelvin- Voigt model for viscoelastic material 
|46| . It consists of a dashpot with viscosity r] and a spring of stiffness k connected in parallel. 
These two modules represent viscous and elastic properties of the material, respectively. 
The Kelvin- Voigt model is the simplest viscoelastic model which in the stationary state is 
determined by elasticity, in contrast to the Maxwell model, which flows in the stationary 
limit. Thus the Kelvin- Voigt model is the appropriate choice for stress fibers, which can carry 
load at constant deformation over a long time. In order to cope with the contractile behavior 
of stress fibers, we introduce an additional contractile element that represents the activity 
of motor proteins. For illustrative reasons we first derive the governing differential equation 
for such a viscoelastic and contractile element {vec- element) before we proceed to the stress 
fiber model. This ansatz is similar to the two-spring model which we have introduced before 



to explain the physical aspects of rigidity sensing on soft elastic substrates [47 1. 

The considered vec-element is depicted in Fig. |5l The properties of the additional con- 
tractile element is given by the specific force- velocity relation of the molecular motor. For 
simplicity, we use the linearized relationship: 

Fmiv) = Fstall (l - -] (5) 



^0, 

Fm is the actual force exerted by a motor moving with velocity v. Vq is the zero-load 
velocity and Fstaii is the stall force of the motor, that is the maximal force allowing motor 
movement. In the following description for the stress fiber, the stall force will depend on 
the active fraction of the myosin, compare Eq. ( |T3l) . such that myosin contractility may vary 
spatially. The sum of all internal forces, F, exerted by the vec-element reads 

= -ku-F^ (6) 

The crucial point is that the motor force Fm{v) is related to the contraction velocity v which 
relates to the displacement u{t) simply by 

du , , 
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The minus sign comes from the fact that the directed motor movement causes a relative 
shding of the anti-parallel orientated actin filaments leading to contraction of the element 
with velocity v. Thus the displacement u{t) becomes negative upon motor activity. If there 
are no external forces, all internal forces have to balance, F = 0, and we arrive at the 
following differential equation for the displacement u(t): 

{v + ^)^ + ku = -F,,^u (8) 
fo dt 

The term originating form the motor activity decomposes into two contributions. First it 
increases the viscosity of the damping term (dissipative signature of the motor) and second 
it contributes a constant force to the inhomogeneous part of the equation (contractility of 
the motor). Assume vanishing initial displacement M(t = 0) = 0, the solution to equation 
Eq. ([H]) is given by: 

u{t) = -Ef^ {I - e-"^) (9) 

Here r = rje/k corresponds to the relaxation time of the Kelvin- Voigt model with an effective 
viscosity rj^ = rj + Fstaii/vQ. Thus the system appears to be more viscous for stronger 
and slower motors and subsequently the motor activity slows down any relaxation due 
to perturbations of the system. The steady state deformation u{t oo) = —Fstaii/k is 
determined by the ratio of stall force and spring stiffness. 



IV. CONTINUUM VERSION OF STRESS FIBER MODEL 

We model a stress fiber as a string of vec-elements whose viscoelastic and contractile 
properties may vary spatially. For example the contractile strength of the vec-elements 
will vary spatially because the biochemical signal will cause different myosin activation 
levels along the fiber. Starting from a discrete description depicted in Fig. [6] we derive the 
governing continuum equation. Similar to the preceding discussion, the force -F„ on the site 
n is the sum of spring forces, viscous drag and the forces built up by the motor proteins: 

d d 

Fn = r]n+l-g^{,Un+l-Un)-r]n-^{Un-Un-l) + kn+l{Un+l-Un)-kn{Un-Un-l)+F,n„+i - -^m„ 

(10) 

Here, we allow that the stiffness of the spring, the viscosity as well as the motor force vary 
spatially. In order to deduce a continuum description we expand the functions r],k,u and Fm 
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at X = na assuming that these functions are smooth within the small distance a. Keeping 
only leading order terms yields: 

Note that the leading differential operator dx acts on 77 and n. The same holds for the second 
term. If k does not vary spatially, it simplifies to kd^. Like for a single vec-element, we argue 
that the motor force depends on the displacement u{x,t). The contraction A„ within 
the n-th element generated by the respective motor is given by A„ = — (n„ — The 
contraction velocity is therefore v{x,t) = A{x,t) ~ —adtdxu{x,t). The found expression for 
the velocity is inserted into the force velocity relation Eq. ([5]) leading to: 

/ add \ 

Fm{x,t) = Fstall{x,t) il + -^^«(^'^) ) (12) 

In contrast to Eq. ([5]), here the stall force is not constant but depends on the phosphorylated 
fraction n{x, t) of MLC along the stress fiber. This comes from the fact that along a myosin 
minifilament and depending on MLC phosphorylation, a larger or smaller fraction of myosin 
heads is able to bind to actin and perform ATP-cycles. The more myosin heads are active 
the larger the maximum force that the bundle can exert to the actin filaments. In our 
model we regard the ensemble of myosins within a cross-section of a stress fiber as one large 
contractile unit with an effective stall force that depends linearly on the active fraction n of 
myosin heads: 

Fstall{x, t) = Fmax n{x, t) (13) 

The effective stall force, Fstaii{x,t), would reach the maximum force F^ax if all myosins 
within this cross section would be working [n = 1). In the following we set Fmax = 50 nN 
which will result in boundary forces exerted by the fiber of about 5 nN which corresponds 
to typical values observed in experiments 45[ . Eq. (fTTil together with Eq. (fT2l) and Eq. (fT3ll 
lead to the final model equation for the stress fiber: 

where the effective viscosity rie{x, t) = ri{x) + Fstaa{x, t)/vo, similar to the findings in Eq. ([HD, 
although here the viscosity varies spatially. Interestingly, only the variation of the motor 
force appears on the right hand side of the equation. As a consequence, a homogeneous 
motor activity will not contribute to the displacements within the string. 
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To obtain some intuition for this equation, assume that the spatially varying stall force 
is given e.g. by the steady state solution of the reaction diffusion system, depicted in 
Fig. m such that Fstaii{,x) = Fmaxnss{x). For this simplifying case where the stall force does 
not vary in time, Eq. (|T4|) can be integrated and the time dependent solution for u{x,t) is 
given by: 



u{x,t) = u{xo,t) + / dx' i dx'u{x' ,tQ)e ^(^') 



Fstall{x' 

ak{x') 



1 — e 



t .If 



+ —e / Fi,it')e<-') dt' ] (15) 

ane^x') Jt, J 

Here, we have set t{x) = rie{x)/k{x), the typical relaxation time with which perturbations 
decay at a certain position. E.g. the initial conditions dxu{x,to) can be regarded as per- 
turbations to the steady state and they decay with exp{—t/T{x)). The three integration 
constants can be identified as the displacement at the left boundary, u{xo,t), the force ex- 
erted to the left boundary, Ff,{t), and the initial strain along the fiber, dxu{x,to). They are 



determined by the boundary and initial conditions. Experiments by Peterson et al. S 



are 



arranged with cells on stiff substrates to which the ends of the stress fiber are connected by 
focal adhesions. Therefore, the appropriate boundary conditions are fixed ends for the fiber, 
namely u{xo, t) = and u{xe, t) = 0. From the second condition one is able to calculate the 
missing integration constant Fh{t) for any initial condition dxu{x,to). The force on the left 
boundary Fb(t) is given as solution to an inhomogeneous Volterra equation of the first kind: 

/ K{t-t')F,{t')dt' = g{t) (16) 



with the kernel 



^« 1 t-t' 



Kit -t')= / dx' —e ^ (17) 

Jxo arie{x') 

and inhomogeneous part g{t) dependent on the initial condition dxU^x^to): 

git) = r dx' f h - e"^l - dM^', to)e~^A (18) 



xo 



ak{x') 



In order to solve the integral equation Eq. (fT6l) for Fb{t), we calculate its time derivative, 
leading to: 

Fb{t) = IItT - t4t ^^(t - t')F,{t')dt' (19) 



K{0) KiO) j,^ 

This equation yields an explicit expression for the initial force at t=0: 

^.(^o) = ll^y ^ (20) 
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By inspection of the kernel Eq. (fT7I) and the inhomogeneous part Eq. (fTSjl one finds that the 
initial force onto the boundary has a finite value, even for the initial condition dxu{x, to) = 0. 
Eq. (fT9l) also yields an iteration rule for the time course of Fb(t), by applying a quadrature, 
where Fb(to) from Eq. ( l20l) is used as a starting value. The solution for Fb{t) is shown in 
Fig. [71 The boundary force is rising from its initial value and then quickly saturates at 
about 4.8 nN. The result for Fh{t) can then be set into the general solution (Eq. (fT5|) ) for 
the displacement u{x,t) along the fiber. Fig. [8] shows this solution, by using the steady 
state activation level for the myosins {Fstaiii^) = Fmax^ssix)) shown in Fig. |4]and assuming 
the initial condition dxu{x,tQ) = as well as the boundary conditions u{0,t) = and 
u{L,t) = 0. Beside the analytical solution, indicated by circles, we also included the direct 
numerical solution of Eq. (HM for comparison. The numerical solution was derived by using 
the MATLAB algorithm "pdepe". The sinusoidal shape of the function u results from 
stronger contractile motors close to the boundaries causing the fiber elements to displace 
into the direction of the boundaries. Hence the displacement u is positive (negative) along 
the right (left) half of the fiber. It is worth noting that the mechanical equilibration of the 
stress fiber occurs within seconds in contrast to the biochemical system which equilibrates 
over minutes. 



V. ROLE OF FEEDBACK 



We already argued that the system of focal adhesions and stress fibers exhibit a closed 
biochemical and mechanical positive feedback loop. Despite this fact the previous results 
were derived under the assumption that the mechanically triggered biochemical signals at 
FAs originate from a constant force. In order to model the full biological system, the varying 
boundary forces have to be fed back into the equation describing the mechanotransduction 
Eq. ([I]). Since the stress fiber model does not include any cross-links (e.g. intermediate 
contacts to the substrate) the tension 7 within the fiber has to be constant and therefore 
equals the boundary forces: 

F6(t) = 7(x = 0) = ar/e(0, t)d,dtu{0, t) + afc(0)a,M(0, t) + ^.^^^^(O, t) (21) 

This relation now connects the biochemical signaling to the mechanical deformation of the 
stress fiber. Thus, the coupled system of reaction equations, Eq. ([1]) to Eq. (jl]), and the 
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mechanical equation Eq. (fT^ have to be solved simultaneously. This can be done numerically 
by using the MATLAB algorithm "pdepe". The whole system of equations and the used 
parameter values are summarized in Tab. [Ill and Tab. IIIIl respectively. 

By doing a steady state analysis we find, that this system of equations exhibit two stable 
steady states for the used parameter values: the first state is characterized by a generally low 
activation level Ugsix) of the myosin motors resulting in marginal boundary forces whereas 
in the second state myosin motors are non-uniformly activated and the exerted forces reach 
a few nN. This bistability is characteristic for a positive feedback system j^. The first "non- 
active" state would correspond to cells that failed to establish mechanical stress whereas the 



second " active" state correspond to cells that are we^ 
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1 adhered to the substrate. In order to 
we start with the system residing in 



simulate the drug experiments by Peterson et al. 
this "active" state and then, at t=0, we perturb the system by turning on the stimulation 
with calyculin. This is modelled by switching instantaneously the inhibition parameter from 
J = 1 to / = 3, thereby omitting the time delay caused by the internalization of the drug. 
The stimulation with calyculin reduces the phosphatase activity and elevates the myosin 
activation level everywhere leading to a quick increase in the boundary forces exerted by 
the stress fiber. The time course of the force exerted onto the focal adhesions is shown 
in Fig. [9l Subsequently, the positive mechanical feedback triggers additional signaling at 
focal adhesions activating myosin motors preferentially at the cell periphery. This results in 
strong spatial gradients in myosin motor activity, see Fig. [TTl The strong peripheral motors 
then contract the fiber to the cost of the central regions where the fiber has to elongate. This 
can be further analyzed by using the numerical solution for the displacement u{x,t). The 
steady states of the displacement u before and after stimulation with calyculin are shown 
in Fig. [TOl The stimulation strongly increases the displacement along the fiber resulting 
in substantial contraction of the fiber close to the boundaries but in expansion around the 
cell center. This finding becomes more apparent in the shown striation pattern calculated 
from the displacement after stimulation (upper string) compared to the striation pattern of 
a completely undistorted fiber (lower string). The bands close to the boundaries have been 
contracted whereas the bands around the center have been expanded, compare also Fig. [131 
We have to stress that the presented stress fiber model is continuous, thus the model cannot 
distinguish between a-actinin bands or MLC bands. The color code in Fig. [TOland Fig. [T3lis 
therefore arbitrary. We also derive the local relative change of density within the fiber which 
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is given in general as the negative trace of the strain tensor Sprei = (p — Po)/Po = —i^iuij). 
Since the model is one dimensional this simplifies to 5prei = —dxu{x,t), plotted in Fig. [T2l 
The local relative change of band width at a certain position within the fiber, is then 
simply given by: {w{x,t) — wq)/wq = —6prei{x,t) = dxu{x,t). The figure shows that the 
inhomogeneous motor activity causes a contraction of the bands up to about 55% close to 
the fiber ends (the relative change in density is positive), whereas the pattern expand up to 
15% at the middle of the fiber (the relative change in density is negative). 

The experimental time course data for the sarcomer length shown in Fig. [2] is intrinsically 
averaged over a certain area in the peripheral and central regions of the cell. In order to 
compare the model results with the experimental finding we therefore define central (center 
±10/im) and peripheral (edges ±10/im) regions of the cell, indicated by vertical lines in 
Fig. [121 The expected sarcomer length at a certain position along the fiber is given by 

w{x,t) = wq + u{x + Wo,t) — u{x,t) or w{x,t) = wq {1 + dxu{x,t)) for wq <^ L (22) 

In the following analysis this measure is averaged over the defined central and peripheral 
regions, respectively. The deduced time courses for the mean pattern bandwidths in the 
distinct regions are shown in Fig. [131 The expected steady state striation patterns are illus- 
trated as insets. Upon stimulation with calyculin, the peripheral mean bandwidth shrinks 
from its initial value of about 0.97 fim down to 0.83 fim, whereas in the central regions, 
the bands elongate from about 1.03 fim up to 1.13 pm. Interestingly the inital mean band- 
width at the center and periphery yet differ in the initial unperturbed steady state of the 
cell (1.03 fim compared to 0.97 fim). This results from the fact that the unperturbed fiber 
already exerts moderate forces on to the focal adhesions which results in slight spatial gra- 
dients in myosin activation. These gradients then sharpen upon stimulation with calyculin, 
see Fig. [TTl as well as Fig. [121 The model results agree qualitatively with the experimental 
findings by Peterson et al. 38|] and the quantitative measurements are within the same order 
of magnitude (compare Fig. [2|). It is worth mentioning that the amplitude of contraction 
or elongation of the fiber scales inversely with the fiber stiffness k: the softer the fiber, the 
stronger the mechanical deformation will be. Thus, a lower k value would simply explain the 
reported higher values for sarcomer contraction of about 30-40%. In our calculation we used 



k = 4:5nN/ fim, a value reported by 48|. The experimentally measured equilibration time 



of the stress fiber upon stimulation is about 20 min (Fig. [2]) which compares to about 3 min 
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for the model results. These quantitave difference originate from two model simplifications. 
First, we lump the focal adhesion associated processes into one equation thereby we shortcut 
the activation of ROCK and neglect prior activation steps of e.g. Rho-Gef or Rho-GTPase. 
Considering these steps would cause an additional time delay. Secondly, the stimulation 
with calyculin happens instantaneously in the model omitting the time delay caused by the 
internalization of the drug. Refining the model and eliminating these simplifications will 
further decrease the differences in equilibration times. 

To highlight the importance of the mechanical feedback we also include the expected 
results for a system neglecting this feedback shown as dashed lines in Fig. [TT] and Fig. [121 
Here the homogeneous induction of the drug cause an almost uniform elevation of myosin 
activation within the cell. Slight differences between cell center and cell periphery would 
persist but rather marginally extend (Fig. fTTl) . In fact stimulation leads also here to amplified 
distortions of the striation pattern. However, the changes in bandwidths are significantly 
smaller compared to the system incorporating the feedback (Fig. [T2!) . Thus the closed 
biochemical and mechanical feedback loop is an essential feature required to describe the 
strong distortions of striation patterns upon homogenous drug induction. 



VI. CONCLUSION 



Here we have presented for the first time a mathematical model for the closed biochemical- 
mechanical feedback loop triggering the upregulation of focal adhesions and stress fibers 
which is typical for cell culture on stiff substrates. In regard to the biochemical part, we 
present for the first time a reaction-diffusion model for Rho-signaling from focal adhesions 
towards stress fibers. Our modelling is based on an extensive review of the literature, which 
provides the list of diffusion and reaction constants summarized in Tab. IIIII In regard to 
the mechanical part, we introduced a new model for stress fibers which takes into account 
the special viscoelastic and contractile properties of the sarcomeric units. For a linear chain 
of many such units, we derived a continuum equation which we solved both analytically and 
numerically. Combining the two model parts resulted in a complete model for the feedback 
loop of interest. We found that this feedback loop leads to bistability between weak and 
strong adhesion and to strong spatial gradients in the deformation pattern. In fact, the 
experimental results can be only explained by incorporating the feedback loop. Our model 



18 



shows how couphng of mechanics and biochemistry can be modeled in generaL In the future, 
it might be used to address also other issues in the context of cell adhesion, including cell 
adhesion to soft elastic substrates or to cyclically stretched substrates. 

In our model, spatial variations in the deformation pattern result from inhomogeneous 
react ion- diffusion fields. This prediction should be experimentally tested in the future. In- 



arger at the periphery than 
exactly as predicted by 



deed it has already been reported that MLC-phosphorylation is 
at the center (both before and after stimulation with calyculin) 
our model. Our model does not account for inhomogeneities resulting from spatial variations 
in mechanical properties (e.g. local accumulation of myosin II, actin or a-actinin), although 
in principle the model is capable to describe these. A detailed experimental analysis in- 
cluding detailed measurements of local actin and myosin accumulation would be needed to 
disentangle the relative contributions of biochemical and mechanical factors to the inhomo- 
geneous deformation pattern. 

An intriguing aspect of our model is the way different stress fibers might cooperate inside 
a living cell. Conceptually it is easy to generalize our model to describe a system in which 
many stress fibers share the signaling input and many focal adhesions share the mechanical 
output. However, it remains a challenge to model also the dynamics of the actomyosin 
system if it is not completely condensed into stress fibers. 
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Abbreviation 


full name 


function 


GDP/GTP 


guanosine diphosphate/ 
guanosine triphosphate 


small molecule without and with 
a third phosphate group, energy 
source of conformational changes 


GEF 


guanine nucleotide 
exchange factor 


activates GTPases by exchanging 
GDP for GTP 


GAP 


GTPase-activating protein 


stimulates GTP-hydrolysis, con- 
verting active GTPases (GTP- 
bound) to their inactive (GDP- 
bound) form 


GDI 


guanine nucleotide 
dissociation inhibitor 


binds to inactive form of GTPases, 
the complex is soluble in the cyto- 
plasm 


MLC 


myosin light chain 


subunit of myosin II molecular mo- 
tors, regulates myosin binding to 
actin filaments 


MLCK 


1 ■ 1 j_ 1 ■ 1 ■ 

myosin light chain kinase 


phosphorylates MLC 


MLCP 


myosin light chain phosphatase 


dephosphorylates MLC 


MLCP-P 


phosphorylated MLCP 


inactive form of MLCP 


MBS 


myosin-binding subunit 


subunit of MLCP whose phosphory- 
lation makes MLCP inactive 


ROCK 


Rho-associated kinase 


ROCK phosphorylates MLCP at its 
myosin-binding subunit (MBS) 


I 


effect of calyculin 


inhibits MLCP from dephosphory- 
lating myosin, thus enhancing con- 
tractility 



TABLE I: Abbreviations and full names of the biochemical components. For each component, a 
short description of its function is given. 
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Model equations 

dROCKjt) _ riFtit){ROCKtot-ROCK{t)) _ V^iROCKjt) , . n 

dt ~ Ki+{ROCKtot-ROCK{t)) K.i+ROCK{t) V^^) 

dMLCP{x,t) _ V.2MLCP-P{x,t) , p. d'^ MLCP{x,t) / r,\ 

dt ~ K.2+MLCP-P{x,t) + ^ dx^ y^"^) 

dMLCP-P{x,t) _ V-2MLCP-P{x,t) , ^ d'^ MLCP-P{x,t) , q\ 

at ~ K^2+MLCP-P{x,t) + y^"^) 

dn{x,t) ^ V3{l-n{x,t)) _ r.zMLCP{x,t)n(x,t) / 

at K3+{l-n{x,t)) K^sl+n{x,t) K'''-'-^) 

ilVeMM + ujx^t) = -llF^Ux.t) (m5) 

Boundary conditions at x=0,L 

aMLCP{x,t) ,V2 ROCK{t)MLCP{x,t) n r)\ 

dx K2+MLCP{x,t) yuCZ) 

aMLCP-P{x,t) _ V2 ROCK{t)MLCP{x,t) /i o\ 

ax TDp K2+MLCP{x,t) {0C6) 

u{x,t) = stiff boundaries (bc5) 



Abbreviations 

Fstall{x,t) = Fmaxn{x,t) 

r]e{x,t) = r]{x) + Fstaii{x,t)/vo 

Fb{t) = ar]eiO, t)d^dtu{0, t) + ak{0)dMO. t) + F^taiiiO, t) 

TABLE II: Summary of model equations. Eqs. (ml-m4) describe successive biocliemical signaling 
events: (ml) focal adhesion associated activation of ROCK; (m2) and (m3) phosphorylation and 
diffusion of MLCP and dephosphorylation and diffusion of MLCP-P; (m4) regulation of the active 
fraction of the myosins, which is identified with the phosphorylated fraction of MLC. Eq. (m5) is 
the mechanical model equation for stress fibers, where u{x, t) is the displacement along the fiber. 
The boundary conditions for the partial differential Eqs. (m2), (m3), (m5) are given by (bc2), 
(bc3), (bc5) respectively. In Eq. (bc2), (bc3), the upper (lower) sign is valid for the left x = 
(right X = L) boundary. For the sake of clarity, we have introduced the listed abbreviations for 
the stall force Fstaii, the effective viscosity r]e and the force exerted onto the boundary Fj,. The 
presented results have been derived with the assumptions that: (I) The diffusion properties of 
the phosphorylated and unphosphorylated form of the phosphatase are the same, hence D = Dp. 
(II) The viscoelastic properties of the stress fiber do not vary in space, therefore k{x) — > k and 

Fmax 

n{x, t)/vQ. 
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Time dependent reaction variables 


Abbreviation 


Meaning 


Used value 


Reference values 


References 


ROCK 


activated form of ROCK 


. . . 5nM 


> \nM 


[41J 


ML CP 


unphosphorylated form of MLCP 


. . . 1.2/jA/ 


1.2 ± 0.3mA/ 


[49] 


MLCP-P 


phosphorylated form of MLCP 


. . . 1.2fiA/ 


1.2 ± 0.3MAf 


[49] 


n 


fraction of active myosin 


... 1 


[MLC-phosl/ [myosinl 




Reaction constants 


MLCK 


myosin light chain kinase 


O.lfiM 


> lOOriA/ 


[50] 


M 


myosin concentration 


30/jM 


25 . . . 30mA/ 


[51] 


Ki 


Michaelis constant 


5nM 


(no value) 




K-i 


Michaelis constant 


4.7nM 


(no value) 




K-2 


Michaelis constant 


0.1/jM 


0.10 ± 0.01/jAf 


[41] 




Michaelis constant 


15/jM 


(no value) 






Michaelis constant 




52.1 ± 7.1mA/ 


[52] 








34.5 ± 2.8mA/ 


[41] 








IS/jA/ 


[53] 








7.7. . .96.0/jA/ 


[54] 








19 . . . 53/jA/ 


[50] 








20/jA/ 


[55] 


/f_3*M 


Michaelis constant 




lOfiA/ 


[56] 




rate constant 


0.3nM/s nN 


(no value) 






maximum velocity 


1.8nA//s 


(no value) 






rate constant 


2.4 1/a 


2.36 ± 0.10 1/s 


[41] 


R2 


maximum velocity 


4.8^m / s 


7'2 • 




V-2 


maximum velocity 


0.1nA//s 


(no value) 




r3*M 


rate constant 


10 1/s 


2.00 ± 0.36 1/s 


[52] 








3.85 ± 0.095 1/s 


[41] 








5.17 1/s 


[53] 








7.37 . . . 171.3 1/s 


[54] 








70 . . . 100 1/s 


[50] 








4.64 1/3 


[55] 


V3 


maximum velocity 


1.0fiA//3 


r3*MLCK 






rate constant 


21 1/s 


21 1/s 


[56] 


D 


diffusion constant of MLCP & MLCP-P 


l4f.Lm^ /s 


10 . . . lOO^m^/s 


[57] 




effect, react, vol. of FAs 


2.0/xm 


(no value) 




I 


inhibition parameter 


1^3 


(no value) 




Parameters of meelianical model 




stall force 


5071 AT 


(no value) 






maximum motor velocity 


l.Ofim 1 s 


~ 0.1 . . . l^m/ s 


[58] 


a 


vec-element spacing 


l.O^m 


1 .0^771 


[38] 


k 


spring stiffness 




45.77iiV/a 


[48] 


V 


viscosity 


4hnN s 1 iira 


^ rk — 45.7nN s/a 


[48,^ 


L 


fiber length 


50/im 


Si 20 . . . 80fi7n 





TABLE III: Model parameters based on literature searcii. We have set the model parameters 
such that they fit into the reported range. The equation for the phosphorylated fraction of MLC 
is normalized to the total myosin concentration denoted by M. In order to make the involved 
reaction constant comparable to the literature values we give K3, K^s, rs and r_3 scaled with 
M. Eq. ([1]) translates mechanical forces into biochemical activation. For this reason the units of 
the rate constant ri are given as [nM/snN]. The typical relaxation time, r, of stress fibers is 
of the order of a few seconds [59| therewith we roughly estimate the viscosity value as 77 ~ rfc 
where we use r = Is. Myosin activation by calyculin is modelled as competitive inhibition of the 
phosphatase. The inhibition parameter / is switched instantaneously from / = 1 to / = 3. Some 
of the reported values have been measured for the interactions of protein fractions and not for the 
native proteins. Furthermore, the experiments have been done on proteins extracted from different 
species. 
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a-actinin v > > > > >t-> actin filament 



FIG. 1: Cells adhere to the extracellular matrix by integrin-mediated contacts called focal adhe- 
sions. These contacts are the anchor points of stress fibers, which are actin filament bundles held 
together by the crosslinker molecule a-actinin and myosin II molecular motors. The myosins are as- 
sembled in myosin mini-filaments. Due to myosin motor activity stress fibers are under tension and 
exert forces to focal adhesions. This mechanical stimulus initiates biochemical signals (Rho-signal) 
that originate from focal adhesions and propagate into the cytoplasm, altering in turn myosin 
activity. Therefore the system of focal adhesions and stress fibers are connected by a biochemical 
and mechanical positive feedback loop (inset). The spatial part of our model is one-dimensional 
with one stress fibers extending between the two focal adhesions at x = and x = L. 
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FIG. 2: (A) Peterson et al. [33] studied the deformation of stress fibers in fibroblasts by using 
fluorescently labelled a-actinin (in green) and myosin light chain (in red). These two components 
arrange sequentially along the stress fibers and thereby form regular striation patterns. Myosin 
contractility was stimulated with the drug calyculin A. Then an inhomogeneous striation pattern 
results: stress fibers contract at the cell periphery (B) but expand at the cell center (C). These 
results were quantified by time course measurements of the mean pattern bandwidths in the re- 
spective regions. Compared with the control (D), stimulation of contractility leads to very strong 
spatial gradients (E) on the time scale of tens of minutes. 
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FIG. 3: Signaling pathway that controls myosin contractility depicted in its appropriate spatial 
context: mechanical cues are transduced to various biochemical signals at focal adhesions, however 
the precise mechanisms have not been resolved yet. One possible mechanism is that a Rho-GEF 
is activated by a mechanosensitive process at focal adhesions. Rho-GEF then promotes Rho-GTP 
loading and subsequent complexation with Rho-associated kinase (ROCK) which gets activated. 
Active ROCK is able to phosphorylate myosin light chain phosphatase (MLCP) at its myosin- 
binding subunit (MBS). MLCP and MLCP-P are freely diffusible in the cytoplasm and thus can 
reach the myosins in the stress fibers. Increased phosphorylation of MLCP to MLCP-P by ROCK 
effectively leads to increased phosphorylation of myosin light chain (MLC), thus increasing myosin 
contractility. 
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FIG. 4: Spatial dependence of the active myosin fraction 77,(3;, t) at four different time points 
t G {30s, 60s, 120s, 180s} as well as for the steady state, nss{x): We solve the biochemical model 
implying the artificial initial conditions that all components are at zero activation level, but set 
the boundary force to 5nN. Because of this mechanical trigger at focal adhesions, MLC gets 
preferentially activated at the boundaries via the Rho-pathway which leads to a steady increase 
of the myosin activation level. Due to diffusible compounds in the Rho-pathway, the increased 
activation level is smoothed out towards the center of the cell. 
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FIG. 5: The viscoelastic and contractile element (vec-element) consists of a spring of stiffness k, a 
contractile module of contraction force and a dashpot of viscosity ij that are all connected in 
parallel. Spring and dashpot taken alone are the usual Kelvin- Voigt model of viscoelastic material. 
The properties of the contractile module are characterized to first order by the linearized force- 
velocity relation. 
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FIG. 6: We model a stress fiber as a string of vec-elements such that the spring stiffness kn, 
the viscosity rjn and the motor force can vary spatially. E.g. the latter will vary spatially 
due to different myosin activity at the periphery compared to the center. The displacement of a 
certain site n is denoted by u„. The force F„ onto this site n consists of elastic, viscous and motor 
contributions, compare Eq. (jlOp . In the following we assume fix boundary conditions, namely 
u{0, t) = and u{L, t) = where L is the total length of the fiber. 
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FIG. 7: Time course of boundary force Fh{t). The solution of the Volterra equation, Eq. (jl6p . 
indicated by circles, was calculated by applying an iteration rule further explained in the main 
text. The dashed line is the boundary force deduced from direct numerical solution of Eq. (jl4p 
which we include for comparison. For the assumed initial condition, dxu{x, 0) = 0, the boundary 
force increases from its non-zero value at t = 0, given by Eq. ()20p and quickly saturates at somehow 
larger values. 
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FIG. 8: We first analyze the solution of the mechanical equation Eq. (jl4p by assuming a steady 
state myosin activation level, nss(x), shown in Fig. HI For this simplifying case where the myosin 
activation level is not time dependent, Eq. ()14p can be solved both numerically and analyt- 
ically. Fig. [8] shows the analytical solution, Eq. (fT5|) . indicated by circles whereas the direct 
numerical solution of Eq. (jl4p is indicated by dotted lines. The solution is given at time points 
t € {0.1s, 0.2s, 0.5s, 1.0s, 3.0s} as well as for the steady state Uss{x), assuming the initial condition: 
dxu{x, 0) = and the boundary conditions: u{0, t) = and u{L, t) = 0. 
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FIG. 9: At i = calyculin is added to the system. Then myosin motors along the filament get 
activated and further increase the tension within the fiber. The time course of the force Ff,{t) 
transduced to the boundaries is shown. 
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FIG. 10: Steady state solution of the displacement u{x) along the fiber before (dashed line) and 
after (solid line) stimulation with calyculin. The stimulation strongly increases the deformation of 
the fiber resulting in substantial distortion of the expected striation pattern (upper line) compared 
to the striation pattern of a completely undistorted fiber (lower line). The upper striation pattern 
was calculated from the displacement data. The bands close to the boundaries have been contracted 
(about 55%) whereas the bands around the center have been expanded (about 15%), compare also 
Fig.ia 
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FIG. 11: Steady state profile of tlie active myosin fraction before stimulation with the drug (dash- 
dotted line), after stimulation including the mechanical feedback (solid line) and after stimulation 
but neglecting the mechanical feedback (dashed line). For the latter case, homogeneous induction 
of the drug cause an almost uniform elevation of myosin activity. Slight differences between center 
and periphery persist but rather marginally extend, whereas the closed feedback system result in 
an amplification of the spatial differences of myosin activation. 




FIG. 12: Relative change in density along the fiber in the steady state before stimulation with the 
drug (dashed-dotted line), after stimulation with the drug including the feedback (solid line) and 
after stimulation but neglecting the mechanical feedback (dashed line). Positive values correspond 
to a compression of the fiber, whereas negative values indicate elongation. In case of the closed 
feedback (solid line) the fiber strongly contracts close to the boundaries up to 55% but elongates 
at the center to about 15%. In order to compare the model results with the experimental findings, 
we arbitrarily define central (center itlO/um) and peripheral (edges itlO/im) regions of the cell, 
indicated by vertical lines. 
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FIG. 13: Time courses of the mean bandwidth of the fiber in the center (upper curve) and in 
the periphery of the cell (lower curve): The shown values are averages of the bandwidths at the 
center itlO/im for the central region and at the edges itlO/xm for the peripheral region. The defined 
intervals are also shown in Fig. [T2j The expected steady state striation patterns for the two distinct 
regions are shown as insets. These results agree qualitatively with the experimental findings by 
Peterson et al., shown in Fig. [2j 



